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Abstract 

►^ ' In this paper we introduce a projective invarinat measure on the 

■^^ I special unitary group. It is directly related to transition probabilities. 

^^ ' It has some interesting connection with convex geometry. Applications 

^~^ \ to approximation of quantum circuits and entanglement are given. 

o 
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The concept of approximation of one operator by another is important in 
many branches of physics and mathematics. For quantum computing and 
Ch ' information in particular, the approximation of unitary or more general op- 

2 ' erators is crucial for implementation of quantum algorithms. In any approx- 

qh| imation theory the notion of distance is essential for a quantitative estimate 

of the accuracy of approximation. In the approximation of linear opera- 
tors on a Banach( or Hilbert) space, the usual distance function or metric 
is induced by the norm. This approach is particularly useful if we restrict 
d I to affine subspaces of the space of operators. This is because the norm on 

the ambient space induces a norm on the space of operators and the metric 
is defined in terms of the latter. However, if we restrict to some subset of 
operators which may not constitute a subspace, the norm induced metric 
may not seem very natural. For two operators A and B the difference A — B 
whose length defines the distance between A and B may take us outside the 
subset. But the concept of a metric does not depend upon algebraic opera- 
tions. In particular, if the relevant subset is a group we are often interested 
in invariant metrics. That is metrics that remain invariant under left (right) 
translations by the group operations. Of course, invariant metrics are known 
to exist for any compact group. 



We may also view the problem of approximation of an operator from 
another perspective. Informally, one could say that a sequence of unitary 
operators C/„ converges to some unitary operator U if, for any given state 

1 a) the expectation values of the sequence Un converges to that of U . In 
fact we will adopt a weaker criteria. Namely, that they converge in certain 
probability measures. This also turns out to be equivalent to convergence 
in the operator norm. But first, a brief synopsis of the paper. 

In Section II, I introduce a metric on the group of unitary operators 
acting on a Hilbert space. The discussion will be confined to an arbitrary 
but fixed Hilbert space, mostly finite dimensional. Many algebraic and ge- 
ometric properties of the metric are proved. Some interesting connections 
with 2-polytopes (polygons) are discussed. I then derive relations with other 
metrics. In particular, equivalence with convergence in the operator norm 
is shown. I also discuss approximation of quantum circuits. 

In Section III, the definitions are extended to the case when the ambient 
space has a tensor product structure. Some connection between quantum 
state entanglement and the convex geometry of the previous section is ex- 
plored. I analyse bipartite entanglement from a different perspective. I also 
discuss extensions to the difficult case of multipartite entanglement. 

I make some concluding remarks about some aspects not covered in the 
paper which will be investigated subsequently. 

2 An invariant metric on the special unitary group 

First, let us fix some notation. In the following, 7i will denote a complex 
Hilbert space with a fixed inner product < , > . I also use Tin to denote a 
space of dimension n. In the following the dimension of all spaces under 
discussion will be assumed to be finite, unless specified otherwise. Let Un 
denote the group of unitary opertaors in B{7in), where the latter denotes the 
algebra of linear operators on Tin- The corresponding subset of hermitian 
operators will be denoted by L{TCn)- The special unitary group SUn C Un is 
the subgroup of opertaors with determinant 1. I use the standard notation 
C and M for the field of real and complex numbers with usual topology. In 
C"" the standard inner product is used. Thus, if a = (xi, . . . , x„)"^ and (3 = 
(yi, . . . ,yn)'^ G C", where A^ denotes the transpose of the matrix A, then 

{a\P) = ^xlVi 

i 

In this section 7i will denote a complex Hilbert space of dimension n. 
The Hilbert space norm induces a norm A ^ ||yl||, on the space of operators 



on Ti, defined by, 



||j4|| = max \\Aih\ = max (<i>|j4| 

U\\=i 1^1,101=1 



If A is normal then ||^|| = max{|A||A an eigenvalue of A . These and other 
properties of the norm |||| may be found in T. The norm induces a metric 
on the space of operators. 

Let ip Gli. he a unit vector. For any pair of unitary operators U,V £ Un 
define 

D^{u,v) = {i-mu^vm^f''' (1) 

Some of the obvious properties of D^ are given in the following theorem. 

Theorem 1 The function D^ satisfies the following for any unitary opera- 
tors U,V. LetW = UW. . 

i. < D^{U,V) < 1. The first equality holds iff ip is an eigenvector of 
W. 

2. D^ is left invariant with respect to group multiplication in Un ■ Thus, 
for any X GUn 

D^{XU,XV) = D^{U,V) 

Furher, D^{UX,VX) = Dx^- 

3. D^ is symmetric. That is, D^{U,V) = D^{V,U). 
4- D^ satisfies the following 

1 



^wiu - e^^vm' < Diiu, V) <\\{u- vm' 

for some real x 



(2) 



Proof: The first assertion follows from the Cauchy-Scwartz inequality and 
the second and third from the definitions. We prove the last one. Note that 

||(C/ - e'''V)i;f = {iPlilU - V)\U - V))iP = 2(1 - Re(^|C/V|V') 

where Iie{z) is the real part of the complex number z. There is a real 
number x such that {iplU^Ve^^lip) is positive. Then 1 — {ip\WVe''^\ip) < 
l-\{tp\UWe'''\i/:)\'^ == D'^{U, V) and the first inequality in|21follows. More- 
over, 

Diiu,v) = (1 + Mu^vima - mu^v\m 

<2(l-Re((V'|f/V|^)) 



The last equality follows from the fact that |(i/'|f7Ty|^^| < i xhe theorem 
is proved. □ 

As a simple illustration of the use of D^ let us look at the Grover search 
algorithm. A class of quantum algorithms which includes the Grover algo- 
rithm may be reasonably described as an approximation of a unitary opera- 
tor by another with high probability. In the search problem we are required 
to find a state labelled by a non-negative integer a say, from an unordered 
collection of such numbers. The corresponding unknown operator we want 
to approximate may be taken to be Ua, the operator that interchanges the 
basis states |0) and \a) leaves the rest unchanged. Then the Grover algo- 
rithm constructs a circuit, represented by a unitary operator Ga , such that 
Z?|o)(C/a, Ga) < e for some error parameter e > 0. 

Definition 1 For any two unitary operators U, V define 

D{U,V) = maxD^(C/,y) = l-min||V'|| = l\{ip\U^V\ilj)\^ (3) 

Call D(U, V) the u-distance between U and V. 

I have used max{ rain) instead of sup{ inf) in the above definitions since 
the unit sphere Sn-i = {ip \ fV'll = 1 is compact and the respective limits 
are attained. Now, define F^{A) = (^|yl|^) for any operator A. Then, 

D{U,V) = l-min\F4U^V)\^ 
\\ip\\=i 

For an operator A, the set ^(^4) = {F^(^)|||-(/'|| = 1} is called the field 
of values or numerical range of A. It is a well-studied concept in linear 
algebra |2j. We therefore have the first geometric characteriztation of the 
■u-distance. Recall that for a metric space (M, p) with metric p and for 
X € M and K C S, the the distance between x and K{ also denoted by p) 
is defined as 

p{x,K) = inf p(x,y) 

If K is compact then there exists y £ K such that p{x, K) = p{x, y). 
1 D{U, V) is the distance of the set F{U'V) from the origin. That is, 
D{U, V) = {1- p2(0, F(C/tl/)))i/2 



Here the metric p used is the standard Euchdean distance in C 

Before hsting the properties of -D(C7, V) let us further investigate its 
geometric meaning. First note that D{U, V) = D{1, U^V). Hence it suffices 
to study the properties of D{1,W) for a unitary operator W. Since W is 
unitary its eigenvalues lie on the unit circle. If Zj = e*'^* is an eigenvalues of 
W , then < Cj < 27r and the angles are read counterclockwise on the unit 
circle. Then writing an aribtrary vector in the basis of eigenvectors of W it 
is easy to see that the numerical range of W is the convex set 

F{W) = {Y,\x^\''e'^^\ ^|x,|2 = l} 

i i 

That is, F{W) is a convex 2-polytope or polygon whose vertices lie on the 
unit circle. I derive below a simple expression for D{1, W). It depends upon 
a elementary geometric result that seems obvious but the proof does not 
appear to be trivial. As I was unable to find a published proof I give an 
elementary detailed one. 

Theorem 2 Let Zi = e*'^*, i = l,...n he the eigenvalues( possibly with 
repititions) of a unitary operator W . Let Ci 's be ordered such that ci < C2 < 
• • • c„ . Let d = max{ | q — Cj | } . Then, 

I IsinfTr)! if the Zj lie inside a semicircle, , , 

D{\,W) = V ^2^' \ : (4) 

I 1 otherwise. 

Proof: The theorem is intuitively obvious. Let C be the the arc connecting 
the Zj's. If there are two arcs connecting all the eigenvalues let C be the 
smaller of the two arcs. If C contains a semicircle then the origin lies inside 
the polygon of F(Pr). Hence, min /)(0,F(I^)) = 0. That is, L:)(\,W) = 1. 
To make it more precise, observe that D{1,W) = D{l,cW), \c\ = 1. Hence, 
we may assume that ci = 0. If C includes a semicircle then it is either lower 
or the upper semicircle. Suppose it is the upper semicircle. Then there 
must be an eigenvalue on the upper semicirle and another on the real axis 
or below it. Then the triangle joining ci and these two eigenvalues contains 
the origin. 

Next, suppose C lies inside a semicircle, say the upper semicircle. If it is 
any other semicircle then rotate C by multiplying with appropriate number 
c, |c| = 1 so we get all the eigenvalues in the upper semicircle. Then, d = c„. 
Again it is intuitively clear that the line joining the points 1 and e*'^" contains 
the point of the polygon that is closest to the origin. To prove it directly we 



have to show that 

n 

mm{|pi +^Pje*'^'p : < pj and "^Pi = 1} = cos^(c„/2) 

1=2 i 

I win follow essentially geometric intuition to prove this. First, let / be the 
line joining the points 1 and z„. Then, it suffices to prove that the line 
segment I' joining the centre to an arbitray point of the polygon intersects I 
at an intenor piont of I'. That is, for any set {pi, • • • ,Pn|pi > and J2iPi = 
1} the equation 

r^Pie'"' = xe^"^ + (1 - x)e^''", (5) 

i 

has a unique solution with < r < 1 and 0<x<l. Asci = the above 
equation is equivalent to the following pair of real equations. 

n 

r{pi + ^PiCOSCj) = X(l - COSC„) +COSCn (6) 

i=2 
n 

(^Pisincj) = (1 -x)sinc„ (7) 



i=2 



Hence 

_ r{pi + Yj1=2 Pi cos Cj) - COS c„ 
(1 - COSC„) 

^ ^ ^ (Er=2 posing) 



smc 



n 



. (Pi + YTi=2 Pi COS Ci) {YTi=2 Pi sin Ci) _ 1 



(1 — COS c„) sin c„ 1 — COS c„ 



rt— 1 



>r{pi sin c„ + ^ Pj(sin (c„ - q) + sin a) + p„ sin c„) = sin c^ 

i=2 

n— 1 

-r[pisinc„ + 2sin (c„/2)(^pjCOS (c„/2 - a) +p„cos (c„/2)) = sinc„ 



n=2 



Since = ci < C2 < • • • < c^ < vr we have cos (c„/2 — q) > cos(c„/2). 
Hence, on the left side of the last equation the expression 



n— 1 



Pi sin c„ + 2 sin (c„/2) (V] cos (c„/2 - a) + p„ cos (c„/2)) 



n=2 



is greater than or equal to 

n 

pism.Cn + 2sin(c„/2)(^piCos(c„/2)) = sinc„ 

n=2 

Consequently, < r < 1. The fact that r > follows from above since sinx 
is nonnegative in the upper semicircle. From the equation [HI it folios that 

< X < 1 theorem is proved. □ 

The next lemma which is very useful for proving important properties of 
D appears as an exercise in j]^. The proof, essentially geometric in nature, 
is not difficult. It is based on the following fact 1 . 

2 Let A and B he normal matrices and let \A — B\ < e. If the disk D(a, p) 
with centre p and radius p in the complex plane contains k eigenvalues of A 
then the disk I){a, p + e) contains at least k eigenvalues of B. 

1 do not prove it here as a more general result is given in the reference quoted 
above. 

Lemma 1 Let U and V he two unitary matrices whose eigenvalues lie on 
a semicircle of the unit circle. Let the eigenvalues {ai}( resp. {bi}) of A( 
resp. B) he labelled counterclockwise. Then, 

max [at — bi\ < \\A — B\\ 
i 

Next, I prove several important properties of D. 

Theorem 3 For any pair of unitary matrices U,V , the function D{U,V) 
satisfies the following. 

1. Projective invariance 

For any complex number c of modulus 1, 

D{U, cV) = D{cU, V) = D{U, V) 

2. Nonnegative 

< D{U, V) < 1 and D{U, V) = iff U = cV, \c\ = 1. D{U, V) = 
'i'ff U = cV for some complex number c with |c| = 1. 

3. Symmetry 

D{U,V) = D{V,U) 



4- Triangle inequality 

For any three unitary matrices U, V, and W 

D{U, V) + D{V, W) < D{U, W) 

5. Invariance. D is invariant under left and right translations in the 
group Un- 

6. D{U, V) = 1 iff there is a unit vector a such that Ua and Va are 
orthogonal. 

Proof: The first three assertions are straightforward consequences of the 
definitions. I prove the triangle inequality. Let X = U^V,Y = V^W and Z = 
U^W. Using the notation in Theorem El let F{X) (resp. F{Y),F{Z)) de- 
note the convex polygon spanned by the eigenvalues of X (resp. Y, Z). Let 
R{X) = p2(0,F(X)) and similarly for F{Y) and F{Z). Since, Z = XY we 
have 

D{U,V)+D{V,W)-D{U,W) = {l-R{X)fl'^+{l-R{Y)fl'^-{l-R{XY)Y'^ 

If both DiU, V) = D{1, X) and D{V, W) = D{1, Y) > 1/2 then there is 
nothing to prove. Hence, we assume D{U,V) < inm(D(V,W),l/2). This 
implies that R{X) > ■^. Further, we may also assume that D[V,W) < 
1. Using the projective invariance property we may further assume that 
the eigenvalues are ordered so that 1 is the first eigenvalue of both X and 
Y corresponding to phase 0. The preceding assumptions imply that the 
eigenvalues e**^' of X are such that q < vr/S. Similarly, if e*"^' are the 
eigenvalues of are the eigenvalues of Y then dj < vr as before we order the 
eigenvalues counterclockwise such that ci = di = 0. Then, from Theorem |21 
it follows that D{1,X) = sinc„/2 and -0(1,^) = sinii„/2. Let e*^» be the 
eigenvalues of XY. Using Lemmas we get 

|giftn _gidn| = 2sm{{hn-dn)/2) < \\Y - XY\\ = \\l - X\\ =2sin(c„/2) 

It follows that hn — dn < Cn- Now there are two cases. 

1. Case 1 /in > vr. Then D{l,XY) = 1. Put /i„ = vr + c. Then, 
< TT — dn < On — c < Cn- Hence, 

D{1,X)+D{1,Y) = smdn/2+smCn/2 > sin(i„/2+sin ((vr - d„)/2) > 1 



2. Case 2 /i„ < vr. Then D{1,XY) = sin/i„/2 and 

D{l,X)+D{l,Y) =smdn/2+smCn/2 > sin(((i„ + c„)/2) > smhn/2 = D{l,XY) 

In the last inequaUty I have assumed that dn + Cn < vr. Otherwise, 
sind„/2 + sinc„/2> 1 > D{1,XY). 

We conclude that the triangle inequality is valid. 

If W is any unitary matrix then invariance with respect to translations 
means 

D{U, V) = D{WU, WV) = D{UW, VW) 

This follows from 

mm{il:\U^V\i)) = mm {ip\WU'^WV\'4j) = mm {'4j\UWWW\i}) 
U\\=i ||t/'||=i ||i/'||=i 

To prove the last item in the theorem observe that 

D{U, y) = 1 iff min (V'|C/V|V') = 
U\\=i 

That is if and only if there is some unit vector a such that (a|C/Ty|Q,^ = q. 
But then Ua and Va are orthogonal. □ 

The triangle inequality for D is proved using the estimates in terms of 
eigenvalues |21 This implies that all the eigenvectors are available so that 
they may form a basis. But if we restrict to some invariant subspace of the 
full space a complete set of eigenvectors may not be available. I therefore 
give an alternative proof of the triangle inequality. 

Lemma 2 Let Un act on some finite dimensional Hilhert space V and let 
D he defined as 

D{U,V)=mMl-mU^Vm^f/^ 
IIV-II 

Then D satisfies the triangle inequality 

D{U, V) < D{U, W) + D{W, V) 

Proof: Let R{U) = min||^||=i |(^|f7|^)p. From the definitions it follows 
that we have to show that 

R{U^W) + R{W^V) - R{U^V) - 2[(1 - R{U^W)){1 - R{W^V))]^/^ < 1 



Since C/V = U'^WW^'V it suffices to show that for any U,V eUn 

R{U) + R{V) - R{UV) - 2[(1 - R{U)){1 - RiV))]^/"^ < 1 (8) 

Let a G y be such that R{UV) = \{a\UV\a)\'^ and let {ai = a, 02, • • • , an} 
be an orthonornial basis. Then, 

R{U) + R{V) - R{UV) < \{a\U\a)\^ + \{a\V\a)\^ - \{a\UV\a)\^ 

= \{a\U\a)\^ + \{a\V\a)\^ - \ Y,{ai\U\a,) {ai\V\ai)\' 

i 

The last line follows from the resolution of identity / = ^^ \ai) (qjI . Consider 
the last term. 

|^(ai|C/|ai)(a,|y|ai)|2> 

i 

{\{a\U\a)\\{a\V\a)\ - \ Y,{o^l\U\a^){ai\U\o^l)\? > 
{\{a\U\a)\\{a\V\a)\ - (^ Kai|C/|a,)|2)V2(| ^ |^„^|f;|„^)|2)i/2)2 

> {\{a\U\a)\\{a\V\a)\ - (1 - \{a\U\a)\^f'\l - \{a\V\a)\^f/^f 

= 1 - R{U) - R{V) + 2R{U)R{V) - 2[R{U)R{V){1 - R{U)){1 - RiV))]^/"^ 

We use Cauchy-Schwartz inequality for getting the fourth line. Using this 
result we get 

R{U) + R{V) - R{UV) < 

2{R{U) + R{V)) - 1 - 2R{U)R{V) + 2[R{U)R{V){1 - R{U)){1 - R{V))]^'^ 
Hence to prove the inequality [H] it suffices to show that 
R{U)+R{V)-R{U)R{V) < l+{{l-R{U)){l-R{V))y/\l-{R{U)R{V)y/^) 

Since < R{U),R{V) < 1 the above inequality follows from (1 - R{U)){1 - 
R{V)) > 0. □ 

The special unitary group Slin may be viewed from two different perspec- 
tives. First, as a subgroup of the unitary group consissting of unitary matri- 
ces of order n with determinant 1. The second point of view is to consider it 
as factor group. Thus, let 2?„ consist of all constant unitary matrices. That 
is matrices of the form e^^In, c real, where Itv is the unit matrix. Then, 
SUn = Un/T>. The algebraic isomorphism is also a topological homeomor- 
phism. We note that D{U, V) is constant on the cosets of T>. Thus we get 
the following corollary. 

10 



Corollary 1 The function D{U, V) : SUn x SUn -^"^ is a metric on SlAn- 

Now consider the tensor product C™ ® C" amd the action of SUn ® SUn 
on it. The eigenvalues of an operator of the form U ®V are given by UiVj, 
where ttj(resp. Vj) are eigenvalues of [/(resp. V). Then, using theorem [21 we 
can show that 

D{Ui®U2,Vi®V2) = 

Tnin{l,D{UuVl)^/l - D^{U2,V2) + D{U2,V2)^/l - D^{Ui,Vi)) 

(9) 

Let us now compare the w-distance defined above with the standard distance 
induced by the sup-norm. Let U,V ^Un ■ Then, 

\\u - v\\ = \\u{i - uW)\\ = \\i - uW\\ 

Write W = U'V. Since I—W is normal ||1— iy|| = max{|A||A an eigenvalue of 1- 
W}. If {e^'^'^lk = 1, . . . ,n} are the eigenvalues of ordered counterclockwise 
so that < ci < C2 • • • < Cn < 27r. Let Ck be the phase angle that is closest 
to vr. Then, 

||[/-l/|| = 2sin(cfc/2) (10) 

Note that, for real x, \\I — e*^/|| = 2| sinx/2|. Thus, if | sinx/2| is relatively 
large then e*^/ can not be close to /. But the operator e*^/ is simply 
multiplies all the states by a constant phase and hence leaves the projective 
space of quantum states invariant. We see that the distance induced my the 
sup-norm does not have projective invariance. For example, the operators 
/ and — / have maximal distance( =2) between them. The same situations 
exists for the distance induced by the Probenius or trace norm on matrices 
defined by |^|_f = 'Tr{A'^A). The u-distance however has manifest projective 
invariance. This is one of the reasons for the introduction of the metric 
D{U,V). However, the latter is closely related to the standard distance. 

Proposition 1 Let U and V be unitary opertors. Then 

D{U,V) = ^\\U-e^^V\\ 

for some real x. 

Proof: We have ||C/ — y|| = ||1 — VF|| where W = U'^V . First, assume 
that the eigenvalues of W lie on a semicircle. By multiplying W with ap- 
propriate factor e*^ we may assume that it is the upper semicircle and that 



11 



1 is an eigenvalue of W. If we order the eigenvalues e^^^ = 1, • 
terclockwise so that = ci < C2 < . . . < vr, we have D{U, V) = sin (c„/2)( 
see Theorem 12) and \\U — e^^V\\ = 2sin(c„/2)( equation I1U|). The lemma 
follows. 

In the second case if the eigenvalues of W span an arc which includes a 
semicircle then D(U, V) attains its maximal value 1 and we multiply V{ and 
W) with a factor e*^ such that -1 is an eigenvalue of W. Then ||C/ — e"'y|| = 2 
and the proof is complete. □ 



3 Examples and Applications 

The metric D has some obvious physical interpretations. If we think of 
U,V as evolution operators. Thus, we write C/(i,to) and V{t,to) to in- 
dicate the time dependance. If \iPq) is the initial state vector then let 
\tpt) = U{t,to)\'ip{tQ)) and {ip'-f.) = V{t,to)\'tp{to)) be the two vectors at time t 
corresponding to the two evolution operators. Then the square of u-distance 
between U and V at time t is maximum of the quantity (1— the transition 
probability between \ipt) and \il^t)), the maximum being taken over all initial 
states. In the case when the Hamiltonian is independent of time we could 
also visualize this as two consecutive operations on the same system. First, 
the evolution operator U{to,t) followed by V\{t,2t — to) 

A quantum circuit is a unitary operator composed of unitary operators 
of order less than or equal to some fixed number k. If the operators Uc and 
Vc represent two such circuits then we say that Vc e— approximates Uc if 
D{Uc, Vc) < e. We may verify such a claim as follows: 

1. Apply inputs in arbitrary state |a) to Vc- 

2. Apply the output of Vc to the output gates of Uc- Do a projective 
measurement with respect to the pair of projection operators Pa = 
|a)(Q| and P^ at the input ofUc- We are actually applying U^ = U^^. 

3. If the estimated probability for the oucome \a) is > 1 — e for all states 
then Vc e— approximates Uc- 

There are some unsatisfactory aspects to the simplistic approach outlined 
above. First, the estimated probability based on observed relative frequen- 
cies is not the actual probability. This can be rectified by giving upper 
bounds to the diffrence between the two. The requirement that the prob- 
abilities be calculated for all states is impossible to satisfy. It is however 
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sufficient to verify tliat D^{Uc,Vc) is sufficiently small for all vectors ip 
in n + 1 independent bases, where n is the dimension of the underlying 
Hilbert space H. Here, "independent bases" means the following. Let, 
B^ = {ai, . . . , a^}, k = 1, . . . , n + 1 be orthonormal bases in H and let 
P^ = \a^){a^\ be the corresponding projection operators onto af. The 
hermitian operators PJ^ have trace 1. Now a general quantum state is den- 
sity operator, that is, a positive definite operator p of trace 1. If / de- 
notes the identity operator then p — I/n is a hermitian operator with trace 
0. We say that the bases B'' are independent if the traceless operators 
P^ — I/n, i = 1, . . . , n — 1 are linearly independent. Then, the latter span 
the space of traceless hermitian operators. Hence, p — I/n can be written 
as a linear combination of the operators Pf^ — I/n, i = 1, . . . ,n — 1. The 
reader may see |H1 or j^ for details. Then, the state p is uniquely deter- 
mined by the transition probabilities 'Yr{pPJ^). Similarly, the state W ■ p is 
determined by the probabilities Tr{W ■ pPf^). Hence, if Tr(yoPj'^) are close to 
TrCW • pPf') then p and U ■ p will be close. In particular, \i p = \'4)){ip\ is a 
pure state then, KV'lPFl^)^ is close to 1 implying that D{Uc,Vc) is "small". 
We can formalize the abive arguments in case of specific bases( e. g. mutually 
unbiased bases) and get an upper bound on D{Uc,Vc). 

3.1 Quantum Search Algorithms 

In this subsection I discuss application of the u-distance to a class of al- 
gorithms known as quantum search algorithms. The name derives from 
the fact that these algorithms can be adapted to the problem of search 
in an unordered database. I give below a generic description of the algo- 
rithm. Let H = C be a Hilbert space of dimension A^. We are given 
a "standard" B basis in H. Write the elements of B as |1), . . . , |A^)}. 
Suppose that we are given a "blackbox" or oracle unitary transformations 
Oa = I — 2\a){a\, a = 1, . . . ,N. We use a sequence of unitary operators 
Uk, . . . ,Ui interleaved with queries to the oracle. Thus, the quantum circuit 
is given by the unitary operator 

Fa,k ^ UkOaUk-lOa ■ ■ ■ U^Oa (11) 

such that the probability of obtaining the result \a) in a measurement in the 
basis B is greater than 1/2. That is, given the initial state V' 

\{a\Fa,km^>^ + c, c>0 (12) 

where c is positive constant that is independent of A^. We assume that the 
probability distribution over the integers J]\f = {1, . . . ,N} is uniform. This 
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implies that the probabihty of the blackbox operator being Oa is equal(= 
1/iV) for all a £ J]\}. Then we may suppose that ij: is the totally symmetric 
state vector. 



N 

Let Ga be the unitary operator that acts on the "plane" T = Span {tp, |a)} 
leaving all vectors perpendicular to T and permutes ip and \a). Then the 
probability specification^] can be written as 

D^{Ga,Fa^k) > - 

The integer A;, which gives the query complexity is also an estimate of the 
circuit size which is related to time complexity. Let us calculate bounds for 
k. Our method of getting these estimates differs from the original one given 
in [3] and illustrates the use of the concepts introduced ealier. 

Let Vfc = C/fc • • • C/i. IjQi (^ = il) ®ip £ H ® H. Define the following 
operators on H ® H hy their action on the basis {\a) ® \h) in H ® H}. 

F^(|a)®|6)) = |a)0F„,,|6) Vl = I®Vl (13) 

P'i\a)0\b))=e'^'^\a)(dGa\b) (14) 

In the above formula A^ are real numbers to be specified. We may visualize 
the operators F'f^ and P' as controlled operation such that if the first "qunit" 
is a then the Fa^k and Ga are respectively applied to the second. It is easy 
to verify that all the operators are unitary. Now, it follows from the basic 
relation 121 that, for some real x, 

^||F^$ - e'-Vm^ < Dl{Fi Vi) < llF^c^ - F.'^i^ (15) 

We use these relations to get lower and upper estimates of D^{Fl^, V^.). First 
we note that 



F',^-Vm' = \\-^{F;,-Vl)Y,\a)^ 



b\Y.(Pa,k-v,m' 



N 

a 

Using a straightforward calculation( see i6j) we get || X]a(-F'a,fc— Vfc)V'|P < 4A;^. 
Hence, 

DliFkX)<^ (16) 
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Next, since D<s>iF'^,Vl) > 1/V2\\F'^<^ - e'W^<^\\ > l/^/2(||(zP' - e'^Vl)<^\\ - 

\\{zP' — F^)<I>||), where z is a complex number of modulus 1. By choosing 

an appropriate z and using the projective invariance of the function Z?$, we 

get \\{zP' - e^Wl)^ > D^{P',V^) and l/V2\\{zP' - F^)$|| < D^{P',Fi;). 

Hence, 

1 

Using the definition of D$ and the operators, we get 



D^{FlX) > ^D^[P'X)-D^{P\Fl) 



z)*(p',y,') = (i-I^E^"'''H^^V'>l 



2x1/2 



- ^ N ' ^ N' 

In the second step we use Cauchy-Schwartz inequality and the fact that 
X^a K'^l^fcV')^ = 1- On the other hand we also have 

o 

We now define e*'*'" = {a\Vk)'4' /\{a\Vkii)\ if |(a|VA:'0)| ^ and 1 otherwise. 
Then, from equation [El it follows that e"*^"(a|Vfc^) = |(a|Vfe^)| > (1/2 + 
c)^/^. Hence, 

D^{P\F'^<{\-cfl^ 
Combining the estimates for D^{P' , V^) and D^{P' , F^) we get 



As we are only interested in asymptotic behaviour, by taking N large enough 
we have (1 - ^)^/2 - 1 + c > c/i/2. Hence, 

Dl{FlX)><?/^. (17) 

Combining the the two bounds for Dq,{Fl^,Vl?) we get k'^/N > c^/4. That 
is, k = 0{\/N). The complexity of the Grover quantum search algorithm is 
0{\/N) and it is the best possible. 
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4 Some Estimates and Generalizations 

In this section I give some estimates of the metric D in special cases. Let us 
estimate the u-distance for some special unitary operators. The CNOT-gate 
C [H] is a unitary opertor on 4-dimensional Hilbert space C, such that 

C|O)0|1)|O)0|1) 
C|1)0|1)|1)0|O) 

-dimensional orthogonal vectors. It 

D{C,U0V)> 1/2 (18) 

where U and V are arbitrary unitary operators in 2-dimensions. In fact we 
can show more. Namely, that the operator C is at ditsance > 1/2 from 
the subgroup of H^ of SU4 generated by the permutation( swap) operators, 
and the product matrices of the form U <SiV and this distance is maximal. 
That is, the CNOT gate is optimal for entanglement of two qubits, a well- 
known result. We may therefore define a measure on the unitary operators 
SU]\[,N = 2" as follows. Let Hn be the subgroup that leaves the set of 
product states invariant. Then for any unitary operator U, 

PE{U,HN) = D{U,HN) = Jnf D{U,V) (19) 

It is conjectured that Hj\f is generated by single qubit operators and permu- 
tations. 

We have defined the distance D on the group of unitary operators via 
their natural representation, i.e. SUn on C", but we could extend it to 
any action of the group. For example, consider the action of SUn the set 
of density operators: C/ • 7 = UjW . But then, we have to be careful in 
checking the triangle inequality. Moreover, if we try a naive extension to 
the infinite-dimensional case we have to deal with convergence issues. I aim 
to deal with these issues in future. 
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